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THE DYNAMIC RESPONSE OF A FLEXIBLE ROTOR BLADE TO 
A CONCENTRATED FORCE MOVING FROM TIP TO ROOT* 

By John F. Ward 
Langley Research Center 

SUMMARY 

An analysis is presented for an approximate solution of the dynamic response of a 
flexible rotor blade to a radially moving force. The results are interpreted specifically 
in terms of helicopter rotor-blade vibration characteristics associated with tip-vortex 
impingement. These results indicate that low-frequency-mode resonant amplification of 
repetitive-impulse loadings may be encountered at moderate to high tip-speed ratios on 
lifting rotor systems. The mode-resonance conditions are encoimtered in sequence with 
increasing tip-speed ratio. The blade natural modes can be forced in resonance at the 
design natural frequencies, which are usually placed at noninteger multiples of rotor 
rotational speed to minimize blade response to harmonic air loads. The results of the 
present analysis suggest that additional attention may have to be given to the selection of 
blade design natural frequencies in order to minimize the response to a repetitive -impulse 
moving force. 


INTRODUCTION 

The impingement of tip vortices, trailing from helicopter rotor blades upon following 
blades, has been identified as a source of dynamic air loads, vibration, noise, and oscil- 
lating stress. This vortex-blade interaction, illustrated in figure 1, has been treated in a 
number of investigations including references 1 to 10. A unique aspect of this type of 
dynamic loading is the radial movement of the trailing vortex along the following blade. 

The pattern of movement on the following blade is systematic and is basically determined 
by the number of blades in the rotor and the rotor tip-speed ratio /i, which is defined as 
the ratio of the forward velocity Vf to the product of the rotor speed and tip 
radius R. 

Emphasis in the past has been on the configurations of low tip-speed ratio (/i < 0.25), 
where the trailing vortex movement is restricted to the outer portions of the following 

*Material herein was presented, in part, in AIAA Paper No. 69-203, "The D3mamic 
Response of a Flexible Rotor Blade to a Tip-Vortex Induced Moving Force," by John F. 
Ward and William J. Snyder, Feb. 1969. 



blade. It is becoming increasingly evident that the problem may be of concern throughout 
the operational envelope for both "pure" and "compound" helicopter lifting rotor systems, 
as indicated in references 11, 12, and 13. At moderate to high tip-speed ratios (/z > 0.25), 
there is a basic change in the nature of the movement of the trailing tip vortex relative to 
the following blade. The vortex movement is no longer restricted to the outer portion of 
the following blade. The point of vortex -blade intersection therefore moves completely 
across the full length of the blade from tip to root at very high relative velocities. This 
full sweep of the blade by a moving loading is the primary reason for concern. As noted 
in reference 14, when a loading moves along a beam at unique critical speeds, resonant 
situations exist which induce increased beam vibrations. 

The objective of the present paper is to provide some insight into the dynamic 
response characteristics specifically associated with the high-velocity radial movement 
of an isolated loading along the blade. The present analysis identifies some of the reson- 
ant conditions which may exist within the operating envelope of pure and compound heli- ' 
copter lifting rotor systems. A basic premise of the present analysis is the presence 
of a strong trailing vortex in close vertical proximity to a following blade. It is 
acknowledged that this condition may not exist for certain operating conditions such as 
zero rotor lift on compound helicopters oi; cases where combinations of rotor lift and 
propulsive force produce rapid movement of the trailing wake away from the rotor tip- 
path plane. 

This analysis does not treat a fundamentally new dynamic problem, since there are 
existing generalized aeroelastic rotor computer programs which include, within detailed 
aerodynamic and d 3 mamic subprograms, all the ingredients of concern in the present 
analysis. Therefore, this simple analysis merely isolates a particular aspect of rotor 
d 3 mamics and emphasizes the blade response to radially moving forces in hopes of gaining 
some additional insight. 

The loading used in this initial analysis is a concentrated unit force. The prediction 
of the actual vortex-induced loading on a following blade depends upon an accurate defini- 
tion of the free wake of the lifting rotor. The detailed consideration of typical vortex- 
induced blade loadings, such as moment couples or locally distributed loadings, are not 
specifically treated in the present analysis. The moving-force analysis is considered to 
be a reasonable first step for the purpose of illustrating basic trends in blade-response 
characteristics, such as resonant conditions, single -impulse response as compared with 
repetitive-impulse response, and associated blade-root bending -moment time histories. 
The analysis illustrates cantilever -blade -response characteristics, and the implications 
of the results are treated in relation to helicopter rotor vibration. The derivation of the 
equations of motion is presented in the appendix. 
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SYMBOLS 


C(x) 

d 

El 

F(x,t) 

Fd(x,t) 

g 

Ki 

ki 

kiR 

M(x,ti) 

m 

Mw 

n 


frequency parameter related to velocity of moving force, radians 

' R ’ 


per second 


damping coefficient at point x on blade, 


kiR\, 


/ 


Ibf 


( newtons \ 
meter/sec i 


ft/sec ^meter/sec^ 
distance from free end of blade to point of load application, ft (meters) 

blade bending stiffness, Ibf-ft^ ^newton-meters^) 

time-varying force applied to blade, Ibf (newtons) 

time -varying damping force applied to blade, — 

ft/sec ymeter/sec 

gravitational acceleration constant, ft/sec^ ^meters/sec^) 
rotating-blade Southwell coefficient 


normal-mode coefficient 


eigenvalues of frequency equation for simple, uniform, nonrotating blade 

time-varying blade bending moment at point x, Ibf-ft (newton-meters) 

mass per unit length of blade, (kilograms^ 

ft \ meter j 

blade deadweight moment at fixed end, , Ibf-ft (newton-meters) 

A 

total number of modes considered 


P(t) time-varying magnitude of concentrated load on blade, Ibf (newtons) 

q(t) time -varying magnitude of nondimensional concentrated load on blade 

qo nondimensional unit force, where applied force P(0) at t = 0 is equal in 

magnitude to blade weight mRg 

R blade radius, ft (meters) 



t time, sec 

tj instant of time for which blade dynamic response is calculated, sec 

V(x,tj^ time -varying shear force at point X on blade, Ibf (newtons) 

Vf forward velocity of aircraft, ft/sec (meters/sec) 

V relative velocity of radially moving force, measured along blade span, 

positive when moving from tip to root, ft/sec (meters/sec) 

Xi(x),Xj(x) modal displacement of ith and jth mode, respectively, at point x on blade 

X radial coordinate along blade, measured from free end, ft (meters) 

y blade lateral displacement, ft (meters) 


«i 


coefficient used in characteristic modal functions 


6(x - d) Dirac delta function 


nondimensional rate of change of magnitude of applied concentrated unit 
force, 

■JoA 


ju, rotor tip-speed ratio, 

assigned value of modal damping coefficient 

T time interval during which force is on blade, sec 

(;o(t) time-varying mode participation factor, ft (meters) 

p blade azimuth angle measured from downwind position in direction of 

rotation, deg 

Ap azimuth increment corresponding to time of vortex movement along blade, deg 

Q, rotor rotational speed, radians/sec 
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<x) blade bending natural frequency, radians/sec 

time -varying slope of blade deflection curve at point x, radians 

Ai^ 

Bi 

Bi ) frequency parameters defined in equations (A38) of appendix 

Ei 


^i,e,l’* • *?^i,e,4 
^i,e,lj* • -7^i,e,4 


coefficients in equations (A34) and (A35) of appendix for ith- 
mode participation factors; these coefficients are defined in 
equation (A37) and are related to an applied force with a 
linear variation in magnitude with time (e ^ 0) 


C. 


i,o,i. 




coefficients in equations (A34) and (A35) of appendix for ith- 
mode participation factors; these coefficients are defined in 
equation (A36) and are associated with an applied force of 
constant magnitude (e = 0) 


Subscripts: 

C chordwise degree of freedom, blade motions in plane of rotor rotation 

d damped 

F flapwise degree of freedom, blade motions perpendicular to plane of rotor 

rotation 



mode number 

R 

rotating 


A prime indicates a derivative with respect to x. 
A dot indicates a derivative with respect to time. 
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TIP-VORTEX MOVEMENT 

The basic concern in the present analysis is the radial movement of the trailing 
vortex from the leading blade of a helicopter along the following blade. The path of a tip 
vortex deposited into the rotor wake by a given blade will describe a systematic helical 
pattern when viewed from above the rotor disk. The instantaneous point of the intersec- 
tion of the following blade and the trailing tip-vortex helix (point A in fig. 1) depends upon 
rotor tip-speed ratio ji and the number of blades. 

An example of the resulting radial movement of the point of blade-vortex intersec- 
tion along the following blade for a three-blade rotor is shown in figure 2. This figure, 
similar to figures presented in reference 1, illustrates typical variations of the intersec- 
tion point as a function of azimuth position xj/ of the following blade for a range of tip- 
speed ratios from 0.2 to 1.0. The point A on the = 0.2 line in figure 2 corresponds 
to the same point in the illustrated example shown in figure 1. The effects of wake con- 
traction, neglected in this presentation, will cause slight changes in detail but will not 
alter the basic pattern variation with increasing tip-speed ratio. At tip-speed ratios 
below approximately 0.25, the tip-vortex path moves onto the following blade at an azimuth 
angle \f/ of about 80®, moves inboard on the blade, and then returns to the tip. As the 
tip-speed ratio is increased above 0.25, a basic change occurs in the pattern of movement. 
The point of intersection again starts at the tip but continues to move inboard along the 
entire blade until it leaves the blade at the root end, instead of returning to the tip. Fur- 
ther increase in tip-speed ratio results in increasing the relative velocity of the move- 
ment of the point of intersection along the blade. 

The situation will change somewhat with the number of blades in a given rotor, as 
shown in figure 3. An increase in the number of blades decreases the interval between 
the leading and the following blade. This decrease, in turn, increases the tip-speed ratio 
at which the vortex path changes from leaving the blade at the tip to leaving the blade at 
the root. 

The vortex-blade intersection patterns for a three-bladed rotor shown in figure 2 
result in the range of radial velocities of vortex movement presented in figure 4. This 
figure presents the variation in the average velocity of radial movement of the vortex- 
blade intersection point along the following blade. Figure 4 is based upon an assumed 
operating tip speed fiR of 650 feet per second (198 meters per second). The curve 
begins at /Lt = 0.25, which corresponds to the tip-speed ratio at which the vortex begins 
to sweep over the entire blade and moves off the following blade at the root end. Typical 
relative velocities are 250 to 600 feet per second (76 to 182 meters per second) for a 
range of tip-speed ratio from 0.25 to 1.0. 
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ANALYTICAL METHOD USED IN DETERMINING BLADE DYNAMIC RESPONSE 


The dynamic response which is amenable to analytical solution applies to the force 
moving along the following blade at constant velocity. This situation is approached as the 
tip-speed ratio increases above 0.25. For tip-speed ratios below 0.25, the intersection 
point of the tip-vortex path and the blade moves along the blade with a varying velocity 
and results in blade-response equations of motion with time-varying coefficients. Some 
insight into these low tip-speed-ratio responses may be gained from the constant velocity 
solutions, but the practical determination of the varying velocity responses requires the 
utilization of numerical integration techniques. 

The problem of determining the dynamic response of a flexible rotor blade to a 
moving force is analogous to the problem of a concentrated load moving across a simple 
beam, which has been treated by a number of investigators (refs. 14, 15, and 16). The 
principal concern in the past has been with severe vibration induced in bridge structures 
by vehicles moving at critical speeds. The analytical method developed herein treats the 
adaptation of the simple beam analysis of reference 14 to a rotating cantilever beam and 
the interpretation of the results in terms of conventional helicopter-rotor-system param- 
eters. The analysis is readily adaptable to hinged rotor blades with the proper substitu- 
tion of hinged blade-mode shapes, frequencies, and damping coefficients in place of those 
of the cantilever blade. 


Assumptions 

The basic assumptions, which result in an approximation of the rotating rotor blade 
dynamic response characteristics, are as follows: 

(1) The rotor blade has uniform mass and stiffness distribution. 

(2) The mode shapes of the rotating flexible blade are reasonably approximated by 
the uncoupled mode shapes for a nonrotating uniform beam. 

(3) The analysis, developed for the flapwise degree of freedom, can also be applied 
to the chordwise degree of freedom by assigning appropriate values of the rotating -blade 
natural frequencies and damping coefficients. 

(4) The periodic variations in aerodynamic damping caused by increasing forward 
speed are neglected. 

(5) The vortex-induced blade loading is considered to be sufficiently localized in 
relation to the blade radius so that it may be treated as a point load of constant or linearly 
varying magnitude, moving at constant velocity along the blade. (This concept is acknowl- 
edged to be an oversimplification of the real situation in light of such factors as vertical 
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displacement of the vortex in the nonrigid wake, time variations in blade lift, vortex 
decay, and complex blade -wake interaction.) 

(6) Only the tip-vortex trailing from a single leading blade is considered. 

(7) The analysis is performed on the basis of an arbitrary "unit force" which is 
nondimensionaUzed by the weight of the blade. 

Equations of Motion 

The equations giving the undamped and damped, forced and free response for the 
cantilever blade arrangement shown in figure 5 are developed in the appendix. In the 
analysis the force is assumed to move from the blade tip to the blade root at a constant 
velocity v and the nondimensional unit force is assumed to vary as a linear function of 
time. The manner in which this unit force q(t) varies with time is illustrated in fig- 
ure 6. The time interval t, indicated in figure 6, represents the time it takes the force 
to move the length of the blade R at constant velocity v. 

Sample Applications 

A few sample applications are presented and discussed in the following section to 
illustrate the type of blade response that might result from tip-vortex impingement at 
moderate to high tip-speed ratios. The approximate analysis is applied to both the flap- 
wise and chordwise degrees of freedom by using the respective rotating natural frequen- 
cies and assuming that the nonrotating normal mode shapes adequately approximate the 
rotating mode shapes in all examples. The basic characteristics for the cantilever blade 
used in the sample cases are Listed in table 1. 

RESULTS AND DISCUSSION 

Resonant Conditions 

Just as is true for the stationary simple beam, the rotating blade experiences 
resonant response characteristics when specific values of the frequency parameter aj 
in equation (A39) of the appendix become equal to a given mode natural frequency 
Since a^ is a function of the velocity of the moving force and, in turn, of tip-speed ratio, 
resonant conditions will be encountered at specific tip-speed ratios. 

The relative velocity v at which resonant conditions exist is determined as 
follows : 

At resonance 


ai = WiR 
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and equation (A28) gives 


Therefore, 



/kiR' 
V = CiR ^ 


-1 


With the relative velocities determined in this manner, the critical tip-speed ratios can 
be obtained from figure 4 for a three -blade rotor with a tip speed f2R of 650 feet per 
second (198 meters per second). 


Applying this approach to the sample cantilever blade described in the previous 
section results in the critical tip- speed ratios shown in figure 7. Resonant conditions, 
implying increased rotor vibration, are encountered in the flapwise modes (indicated in 
fig. 7 by the circular symbols) at a tip-speed ratio of approximately 0.5. The first three 
flapwise bending modes are encountered almost simultaneously, and the fourth and higher 
modes are encountered in sequence at higher tip-speed ratios. Applying the same 
approach to the chordwise bending degree of freedom results in a predicted resonance 
for the first mode at a tip-speed ratio in the vicinity of 0.7 (shown by the square symbol 
in fig. 7). This example indicates the general trend in resonant points. The points shift 
somewhat with changes in blade frequencies, tip speed, and number of blades. For 
example, most of the hinged-blade mode resonant responses are encountered over the 
same range of tip-speed ratios as shown in figure 7. However, the hinged-blade chord- 
wise rigid-body mode, with the natural frequency a subharmonic of rotor rotational speed, 
can have a critical tip-speed ratio as low as 0.3. 


A decrease in the number of blades at a given tip-speed ratio increases the radial 
velocity of the tip vortex induced moving force. This trend can be observed in figure 3 
by noting the increase in the slope of the vortex-movement curves with decreasing num- 
ber of blades for a constant tip-speed ratio ju. = 1.0. This result means that a given 
resonant radial velocity will be achieved at a lower tip-speed ratio when the number of 
blades is decreased. Conversely, the addition of blades can decrease the radial velocities 
and delay resonance to higher tip-speed ratios. 


The moving force is on the blade during only a portion of each rotor revolution. 
Therefore, the resonant blade motion will not expand continuously as a function of time. 
While the moving force is on the blade, during the time interval t, the equation for the 
forced response mode participation factor at resonance is as follows: 


For aj^ = 0 < tj < t, = 0, and e = 0, 
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(p. 


iR' 


('i)“ 


Qog 


sinh coijjti + Q!i(cosh - cos 


" tl"iR(<^os coii^ti - Q-i sin Wi^tj) 


This equation was derived from equation (A39), which is given in the appendix. Since 
equation (A39) is indeterminate at resonance, ai = this equation was obtained by 
taking the 


lim 

at - ^iR 


f'(%) 

F'(ai) 


where f'(ai) and F'(ai) are the derivatives of the numerator and denominator, 
respectively, of equation (A39). 


Blade -Root Bending Moments 

The following examples and discussion deal with the basic cantilever blade 
described in table 1. Examples of single- and repetitive-impulse blade-root bending- 
moment response are given. Both the flapwise and chordwise degrees of freedom are 
discussed as required for illustrative purposes. The damping used in various cases is 
assigned to cover a range of values considered to be representative of cantilever blade 
designs. These constant damping coefficients are assumed to include both structural and 
aerodynamic damping. The periodically varying damping coefficients which result from 
the effects of increasing tip-speed ratio cannot be treated in this approximate analysis. 

The specific examples considered in the following discussion are listed in table 2. 
The fourth flapwise bending mode, identified in figure 7, is not included in these sample 
cases. 

Single imp ulse.- If the damping is sufficient to reduce the response to insignificant 
levels within one rotor revolution, a single -impulse analysis can be used. This type of 
analysis generally applies in the flapwise degree of freedom where aerodynamic damping 
is significant. Examples of the blade-root flapwise-bending-moment variation with time 
for a single-impulse loading are shown in figure 8, which presents three examples of the 
damped response to the impingement of a moving unit force. These three examples 
represent situations where the concentrated force is moving along the blade at relative 
velocities of 248 feet per second (75 meters per second) when A\p = 150°, 372 feet per 
second (113 meters per second) when A\p = 100°, and 466 feet per second (142 meters 
per second) when Ai// = 80°. These results are considered to be analogous to tip vor- 
tices encountered at tip-speed ratios of approximately 0.3, 0.5, and 0.7, respectively. 
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The single -impulse results shown in figure 9 are for a tip-speed ratio of 0,7 with 
reduced modal damping. The transient response persists for more than 2 revolutions 
and the difference between the two response signatures is the result of using a moving 
force of constant magnitude in one case (e =0), shown by the solid line, and a moving 
force decreasing linearly to zero (e = -1) in time t, shown by the dashed line. 

Repetitive impulse .- The rotating blade experiences a tip-vortex-induced impulse 
loading each revolution. If the combined aerodynamic and structural damping is low, 
the free vibration will persist over a number of rotor revolutions. This response may 
be approximated by the superposition of a series of single -impulse response time his- 
tories, The summation is continued for a sufficient number of rotor revolutions to 
achieve steady-state conditions. An example of low damping can be considered in either 
the flapwise or chordwise degree of freedom. 

Chordwise degree of freedom: The chordwise case is considered here first for 
illustrative purposes. The example shown in figure 10 represents a typical blade-root 
chordwise-bending-moment time history at a tip-speed ratio of approximately 0.3. The 
low level of damping assigned to this response is assumed to represent a typical 
chordwise-mode damping which consists of structural damping and little or no aerody- 
namic damping. The second and higher chordwise modes are not considered because of 
the wide frequency separation between the fundamental and second chordwise modes. 

The solid curve in figure 10 represents the blade-root chordwise-bending-moment 
response to the initial impulse in the first rotor revolution. Because of the low damping, 
the transient response persists over many revolutions. Eventually, the repetitive- 
impulse response will stabilize on a steady- state response represented by the dashed 
curve in figure 10. This steady-state response, after 10 rotor revolutions, appears as a 
distorted once-per-rotor revolution periodic oscillation. 

Increasing the tip-speed ratio for the example shown in figure 10 from 0.3 to 0.7 
results in the root-bending-moment signatures shown in figure 11. Here again the 
chordwise-moment time history is shown for the initial revolution and the steady-state 
case in the tenth revolution. The result shown in figure 11 has a peak-to-peak response 
which is lower than the response for a tip-speed ratio of 0.3 in figure 10. This result is 
true even though the result at a tip-speed ratio of 0.7 is close to a critical tip-speed ratio 
for the chordwise mode indicated in figure 7. This decrease in response results from a 
specific combination of forced and transient response of the blade when rotated at a non- 
integer multiple of the chordwise fundamental bending-mode natural frequency ^in this 
case = 0.7). In addition, the time interval t during which the force is on the 

blade decreases with increasing tip-speed ratio. 

Flapwise degree of freedom: An example of the type of response that might be 
encountered in the flapwise degree of freedom, with low damping in each mode, is shown 
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in figure 12. Both the single -impulse and repetitive -impulse signatures are shown. The 
repetitive -impulse loading again results in a steady- state signature (repeating every 
rotor revolution) that is a mbcture of the forced and transient response of the blade flap- 
wise modes. This combination gives a rather complex response, which would prove dif- 
ficult to interpret realistically on the basis of applying conventional harmonic analysis 
techniques. 


Further Comments 

Generalizations, similar to those for a simple harmonic forced response, are diffi- 
cult to state for the repetitive impulse response. Calculated time histories of blade 
response for a wide range of conditions suggest that the maximum response for conven- 
tional blade designs will occur in the tip-speed-ratio range from 0.25 to 0.35. The ana- 
l 3 dical results also confirm that, as expected, when the blade natural frequencies are 
placed at integer multiples of rotor rotational speed, and the low frequency modes are 
lightly damped, the response to repetitive-impulse loading can build to very large 
amplitude. 

Proper tuning of the rotor -blade frequencies to specific noninteger multiples of 
rotor rotational speed may make it possible to attenuate the blade response to tip -vortex 
repetitive-impulse loadings. This situation suggests the need for some additional atten- 
tion in the selection of design natural frequencies of the blade to avoid the amplification 
of harmonic airloads. In this connection, reference 17 treats some of the fundamental 
aspects of tuning dynamic systems to achieve minimum response to repetitive impulse. 
Also, reference 18 treats a case of helicopter rotor dynamic stability associated with 
response to repetitive impulse. 

The examples presented in this paper are for a cantilever, or hingeless, rotor blade. 
The method may be applied to hinged blades as well. The resulting critical tip-speed 
ratios will be similar to those of the cantilever blade and will be directly related to the 
particular natural frequencies, tip speed, and number of blades of the hinged rotor system. 

A further extension in the analysis would utilize rotating mode shapes in order to 
provide increased accuracy in the response calculations. Each rotating mode may be 
included in the analysis as an expansion in normal modes of a nonrotating blade. The 
expansion method is discussed in reference 19. Some insight into the case of nonconstant 
velocity of the moving force, which applies to low tip-speed-ratio (/i < 0.25) conditions, 
can be obtained by considering classical dynamic response during an acceleration or 
deceleration through critical shaft speeds, as discussed in references 16 and 20. 
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CONCLUDING REMARKS 


An analytical study of rotor -blade djmamic response to impulse loading indicates 
that low-frequency-mode resonant amplification can be induced by the force moving at 
high velocity from the tip to the root of the blade. The practical implication of this phe- 
nomenon is that the situation is analogous in many respects to the impingement of a 
trailing tip vortex upon following rotor blades at moderate and high forward speeds. 
Interpreting the results in relation to tip-vortex movement along a rotor blade indicates 
that blade mode resonance may be encountered at specific tip-speed ratios. These mode 
resonant points are characterized by increased blade vibration and are encountered in 
sequence at increasing tip-speed ratios above a value of approximately 0.25. The blade 
natural modes can be forced in resonance at the design natural frequencies, which are 
usually placed at noninteger multiples of rotor rotational speed to minimize response to 
harmonic air loads. The results of the present analysis suggest the need for additional 
attention in the selection of blade design natural frequencies in order to minimize the 
response to a repetitive -impulse moving force. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 15, 1969, 

721-02-10-05-23. 


APPENDIX 


EQUATIONS OF MOTION 

The cantilever -beam arrangement shown in figure 5 rotates at constant angular 
velocity S2. The x,y-coordinate system rotates with the beam with the origin at the point 
of the undeflected, free end of the cantilever beam. This selection of an unconventional 
coordinate system resulted from the primary concern for studying beam response for 
loadings moving onto the beam at the free end (x = 0) at the initial time (t = 0) and moving 
radially at constant velocity v to the fixed end of the beam at x = R and t = t. 

The analysis proceeds by first developing the equations of motion for a nonrotating 
beam, fi = 0, by the methods of reference 14, and then using these results to provide an 
approximate analysis for a rotating blade. The approximation is achieved by assuming 
that mode shapes for a nonrotating beam adequately represent the mode shapes for a 
rotating blade and that the natural frequencies of the rotating blade can be calculated by 
accounting for centrifugal stiffening effects by the method of reference 17. 

Nonrotating Beam 

The governing differential equation of forced motion of the nonrotating beam is 

my(x,t) + |ily”(x,t[| = F(x,t) (Al) 

The external forces acting on the beam are the viscous damping force and a concen- 
trated force P at a point x = d: 

F(x,t) = F(j(x,t) + Px=d(t) (A2) 

where the viscous damping force is a function of local beam velocity 

Fd(x,t) = -C(x)y(x,t) (A3) 

where C(x) is the local damping coefficient at point x. Therefore, equation (Al) 
becomes 

my(x,t) + C(x)y(x,t) + [EIy"(x,t)] = 

Using the Galerkin method permits the deflection y(x,t) to be represented by a 
superposition of modal deflection fvinctions as follows: 
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APPENDIX 


n 

y(x,t) = ^ <p.(t)Xj(x) (A5) 

3=1 

where the functions Xj(x) are the normal mode shapes satisfying the equation 


|eIXj'’(x)J" - mo)j2x.(x) = 0 


(A6) 


The functions cp^ (t) are the mode participation factors and coj is the frequency of 
oscillation of the jth mode. 

Substituting equation (A5) into equation (A4) gives 

n ^ ^ 

Z [? ^3 " ”^^3 *^3 ^ *^3 ^*^^3 ^x=d(t) 

j=iv. J 


Applying equation (A6) to equation (A7) gives 


11 

I - 


m<; 0 j(t) 


+ C(x)(^:i{t) + mWj2<;(,.(t) 


Xj(x) = Px=dW 


3=1 ^ 

Multiplying equation (A8) by Xj[(x) and integrating over the length of the beam, the 
orthogonality of the modal functions 


(A8) 


\ X;(x)Xi(x) dx = 0 
Jq •* 

permits equation (A9) to be simplified to 


(i ^ 3) (A9) 


n 


mijpiCt) r Xi2(x) dx + y ^i(t) r C(x)Xi(x)X^(x) dx + mw^2(;^7i(t) C Xi2(x) dx 

i^l ^ 

pR. 

= \ Px=d(t)Xi(x) dx 
•^0 


(AlO) 


Neglecting the intermodal coupling due to damping gives i uncoupled differential equa- 
tions for the mode participation factors; 
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^i(t) + 


C(x)Xi2(x) dx 


m Xi^(x) c 
^0 


<pj(t) + = 


f Px=d(t)^iW ^ 


(All) 


m \ A^^(x) dx 

^0 


pR r*R 

The term j Pjj^jj(t)Xi(x) dx can be represented by P(t)6(x - d)Xi(x) dx where the 


function 6(x - d) is the Dirac delta function and 


pR 

\ P(t) 6(x - d)Xi(x) dx = P(t)Xi(d) 

Jq 


Equation (All) becomes 


^i(t) + 


C(x)Xi2(x) dx 

£0 

m r X;^(x) dx 
Jfi 


• 2 P(t)Xi(d) 


m 


Xi2(x) 

Jn 


(A12) 


dx 


The normal functions and frequency equation for a free-clamped beam are developed 
as follows: 

The general solution for the free lateral vibration of a simple uniform beam is of 
the form 


X^ = Cl sin k^x + C2 cos kix + C3 sinh kix + C4 cosh k^x 


For the cantilever beam of figure 5, the following boundary conditions apply: 
Clamped end, x = R: 


Free end, x = 0: 


Deflection 

Xi = 0 

£Ei = o 

dx 

Slope 


d2Xi 

Moment 

1 = 0 


dx^ 


d^Xi ^ 

Shear 

dx3 


(A13) 
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Applying these boundary conditions to equation (A13) and its derivatives and setting the 
resulting determinant equal to zero gives the following frequency equation: 


cos kjR cosh kj^R = -1 (A14) 

This frequency equation has an infinite set of solutions, the first four of which are: 


kjR = 1.876 
k 2 R = 4.694 
kgR = 7.854 
k 4 R = 10.995 


The natural frequencies of the nonrotating uniform beam are given by 


a>i = 



(A15) 


Solving for the constants Cj, . . . ,€4 in equation (A13) gives the following general equa- 
tion for the mode shapes of the nonrotating cantilever beam: 


Xi(x) = sin kjx + sinh kjx + or^^cos kjx + cosh kjx^ 


where 


Q!i = - 


sinh kjR + sin kjR^ 
cosh k^R + cos kjRy 


The first four ai coefficients are 


O'! = -1.36222 
Q !2 = -0.98187 
Qfg = -1.00078 
0:4 = -0.99997 


(A16) 


(A17) 
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The following integral property of the normal modes is utilized to simplify equa- 
tion (A12): 


First, 


Therefore, 


Xi2(x) dx = Q!i2R 

>^0 


<Pi(t) + 


r C(x)Xi2(x) dx 






(A18) 


(A19) 


As in reference 14, if zero initial conditions are assumed, the general solution of equa- 
tion (A19) is 


<;^idW = — f f q(t)Xi(d)e *)sin - t) dt 


“i ^id 0 


(A20) 


which is the superposition, or Duhamel, integral where 


q(t) = ^ 

' ' mRg 


(A21) 


which is the ratio of the magnitude of the applied concentrated force to the weight of the 
uniform beam. The damping ratio in equation (A20) is 


^i = 


C C(x)Xi^(x) dx 

Jn 

2<x)jma!j2R 


(A22) 


Also, 


CO 


id 


(i 2^^/^ 

= (1 - ^i 


(A23) 


and 


Xi(d) = sin(i^jd + sinh^^^jd + 




(A24) 
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If the force P(t) is moving along the length of the beam, it produces vibrations 
that can be calculated by letting v denote the constant velocity at which the force moves 
along the beam, starting at the free end of the beam, x = 0, at time t = 0. At any 
time t, the position of the force on the beam will be 


d = vt 


(A25) 


Under these conditions, equation (A20) becomes 


•^id^*) = \ ^ q(t)Xi(vt)e *^sin - t) 

«r^id *^0 


dt 


(A26) 


where the modal displacement of the point of application of the moving force is 


Xi(vt) = sin ajt + sinh a^t + 0 !j^(cos a^t + cosh ajt^ 


(A27) 


where 



(A28) 


The moving force is considered to vary linearly in magnitude as illustrated in figure 6 
where 


q(t) = qo(l + 1 1 ) 
q(t) = 0 


(0 < t < t) 
(t < t) 


and e is the nondimensional rate of change of the applied force with time; thus, 

^-1 


e = 


Aq(t) 


qo 


(A29) 


The time interval t corresponds to the time required for the concentrated force to 
move from the free end of the beam to the fixed end; namely. 


-I 


(A30) 
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The final equations for the mode participation factors for the damped response to a 
moving force are as follows: 

For 0 < t < T, 


"" a 2*^ ^0 <^id(*l - *) dt (A31) 


4 ‘"id 


For T < t, 


qng r^/ . -*^l‘‘^i(tl-t) 

= - 2 ^- V T <^id(tl - t) dt 


4 ^id 


(A32) 


Rotating Blade 

In order to obtain a first approximation of the response of the rotating helicopter 
blade to a radially moving force, the preceding analysis is used by substituting rotating 
blade natural frequencies. These frequencies are obtained from reference 19, where the 
rotating blade natural frequencies are determined from 


WiR^ = C0j2 + 


The appropriate Southwell coefficients Kj are selected from reference 19 for the flap- 
wise and chordwise bending modes. Also, 


‘‘^idR = 


(' 



^iR 


(A33) 


where the values of effective damping ratio are assigned. 


Using the rotating-blade natural frequencies 
recombining terms from equations (A31) and (A32) 
equations for the rotating blade mode participation 


and expanding, integrating, and 
results in the following approximate 
factors 
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For 0 < tj < T, 


<-'’idR(‘l) ' 


gpe 


1 +Ttl)(Ci,o,l - Ci.o.l) - ?(Ci,e,l - ajH + + 7 ti)(Ci_o ,2 - ^1,0,2) -?(<^i,e ,2 ‘ Ci,e,2^<^os ^iti 




£ .\pT e 7T "^'i^^iR*! Iri 

K ^ 7 *lj^i,0,4 “ 7 ^i,€ 

’T ■ L 

+ 7 tlj^i,0,4 ■ 


■ ■ ^i,e,3 + ^i,6,3^®“^ ^idR*! " [^i,0,2 “ ^i,0,2 + ^i,0,4 " ^i,0,4 

- ?(^i,e,2 “ Cj g 2 + ^i,e,4 ' 4^003 WidRt-\ 


For T < tj, 






2^i^^’idR 


,-^i^iR(tl-^) 


(1 + f 


)Ci.0.1 - 7 ^^sinj^iT - u^idR(ti “ + f)Ci,0,2 - T ^ 2j ^’os|aiT - 


- [^1 ^ f )Ci,0,l - f smjaiT . -'idR(ti - r)J - |(1 + r)Ci q,2 " J , 3 ] " -'idR(tl “ O] 

+ ^ f )C,^o.3 - T C, - jil - f)Ci 0,3 - ? ‘""'^sin ^R|R(tl " ^ 

" <|(1 ^ f )Ci,0,4 - 7 ^i.e , 4 ]^^' - [^(1 + f )Ci.o,4 - T ^i.e ^’idR(tl “ ‘}j 

" ^ ^ ^<^1.0,1 ^ ^1.0,1 ■ ^1.0,3 ^ '^1,0,3 ■ ,1 ^ ,1 ■ *^i.f ,3 + •^’icl 

" 1^1, 0.2 ■ ^i,0.2 ^ ^i,0,4 ~ q, 4 " ^^i,f ,2 " ^i,e ,2 ^ ,4 " ^ 4^Jc 


-^’idRU 


cos ^’idRti 


The coefficients in equations (A34) and (A35) are as follows: 


r. ~ <^iAi 

- Bi - aiAj 

‘-“’1 ■ Ai2 8(2 

'-1,0,1 - ^ 2 ^0 

+ Bi'^ 

^ Aj + o'iBi 

- Aj + cviBi 

^i,0,2 " A 2 n 2 
' Ai'^ + Bi^ 

'-'i,0,2 .2 2 

’ Ai^ + Bi^ 

^‘. 0.3 = 

Et^ + 0j2 

p (1 - «l)Ei 

Ei^ +Di-^ 

p (1 + «0Di 

Ei2 

(1 - «i)Di 

■ Ei2Ti5;^ 


(A34) 


(A35) 


(A36) 
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_ 2Aj^Bi - 

(Aj2 - 6^2) + 

(Ai2 + Bi2)2 


Ci,e,3 - 


_ (Ej2 - D^2 )(i + a.) 


(Ei2n.Di2)2 


Ci,e,4 - 


2Ej^Di(l + Q!i) 


(Ei2+Di2)2 


_ 2A^Bi - (Aj^2 - B^2 ^q,.^ 

J ^ / 9 — 9\2 

(Ai2 + Bi2) 




^i,e,2 - 


_ (Ai2 - Bi2) + 2o!iAiBi 


(Ai2 + Bi2)2 


_ (e^ 2 - p^2)(i _ 


i,e,3 


(Ei2+Di2)2 


2EPi(l - Q>i) 

^’^’^"(Ei2.Di2)2 


) (A37) 


where 


Ai = 


Bi = ai + coidR 
Bi = Oi - 

Ei = ai + 

Ei = -ai + 

= “idR 


(A38) 




Some of the response characteristics can be deduced directly from the undamped- 
blade equations of motion. This equation is given as follows for a constant force: 

For Vi = 0, e = 0, and 0 < tj < t, 


^iR(^l) “ 


Qog 


v(ai2 - WiR2)sinh aitj - (ai2 + WiR2)sin aitj 


Q!i2(ai4 - a>iR4) 

- coii^2^cosh aiti - (aj 2 + coij^2jcos a^ti 


+ a 


+ 2 aiC 0 iR sin c^iRti + 2aioo]ji cos wiRt^^ 


(A39) 
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For T < tj, 


<i^iR(tl) ” „ ' U A' - "iR(tl - ^)] + coslaiT - o)iR(ti - t)^ 

2«i ^iR\^i ■'■ ^iR 1 -’) 


- ■ '")] “i cosj^iT + coiR(ti - 

a-2 0 !i)e^i^ + (l - o!i)e sin WiR(ti - t) 

+ (l + Q'i)e^i'^ - (l - o;i)e‘^^^coiR cos WiR(ti - 
+ 4 (ai^iR sin + o-iWiR^ cos WiRti^ 


(A40) 


ar - ‘^iR' 


If the velocity v of the moving force on a nonrotating beam is very small, the 
equations of motion for t = 0 reduce to the static displacement of the free end of a canti- 
lever beam subjected to a concentrated load at the tip as follows: 


y = 2 ^i(ti)Xi(x) 
i=l 

and equation (A39) becomes, for the forced, nonrotating blade, 

<Pi(tl) = ^ ^ ‘^i^)sin a^tj 

+ Q!i j^aj^ - Wj^^^cosh ajtj - (a^ + coj^^^cos a^t^ 


For tj = 0, V = 0, 
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n 


y=^ <?^i(0)Xi(0) 


i=l 


n 

- V 

2P /^i( 

n 

4Pr3 

y 1 

El 

1=1 (“i") 

_ 4Pr3 

1 

El 

(1.875)4 

= ^^^^(0.08071 - 
El 

= 0.3321 

Pr3 

El 


1 1 

+ - + . 


.) 


The exact value for the static tip deflection is 

Pr3 


y = 


3EI 


With three terms of the series the error is less than 0.5 percent. 

The following general equations give the approximate values of beam deflection, 
slope, bending moment, and shear: 

n 

y(x>ti) « ^ (/>i(ti)x^(x) 

i=l 
n 


0(x,ti) == ^ cPi(ti)x’(x) 
i=l 

n 

M(x,ti) « -El ^ ^i(ti)Xi”(x) 
i=l 
n 

V(x,ti)«-EI^ <Pi(ti)Xi'”(x) 


(A41) 


i=l 
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where 


Xj^'(x) = kjjcosh k|X + cos kiX + o-j^sinh kjX - sin kjx^ 
X’’(x) = ki^jsinh kj^x - sin kjx + Q!j^(cosh kjx - cos kjx^ 
X"’(x) = kj^jcosh kj^x - cos k^x + aj^^sinh kjX + sin kjxJ 


(A42) 
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TABLE 1.- CANTILEVER-BLADE CHARACTERISTICS 


Blade radius, R 17.5 ft (5,35 meters) 

Mass per unit length, m 0.130 slug/ft (6.23 kg/meter) 

Flapwise structural stiffness, (EI)p . . . 4.86 x 10^ lbf-ft2 ( 2 . 01 x 10^ newton-meters^j 
Chordwise structural stiffness, (EI)q . . 2.78 X 10® Ibf-ft^ (l.l5 x 10® newton-meters2) 
Number of blades 3 


Rotor rotational speed 

Flapwise mode natural frequencies, rotating: 

First bending mode 

Second bending mode 

Third bending mode 

Fourth bending mode 

Chordwise natural frequencies, rotating: 

First bending mode 

Second bending mode 


37 radians/sec 

41 radians /sec 
102 radians/sec 
194 radians/sec 
293 radians/sec 

52 radians /sec 
304 radians/sec 


All damping coefficients, fj, are assigned as noted in specific examples 
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TABLE 2.- ILLUSTRATIVE EXAMPLES 


Blade mode 




'Q, 




A^, 

deg 


Figure 


Remarks 


Flapwise bending 


Flapwise bending 



1.11 

2.76 

5.24 




150 

100 

80 



1.11 0.25 1 

2.76 .10 > 0, -1 0.7 80 

5.24 .14 


8 Single impulse 


9 Single impulse 


First rotor revo- 

Chordwise bending 1 1 , 1.41 0.06 0 0.3 150 10 lution and tenth 

rotor revolution 

First rotor revo- 

Chordwise bending 1 1 1.41 0.06 0 0.7 80 11 lution and tenth 

rotor revolution 


Single impulse and 
repetitive impulse 



1 

1 

1 

1.11 

0.13 

1 




Flapwise bending 

3 ; 

2 

1 

2.76 

.05 

) 0 

0.3 

150 

12 


! 

Is 

5.24 

.07 

j 






Point A - Intersection of trailing tip vortex 



Figure 1.- Intersection of following blade with tip vortex trailed from a leading blade at tip-speed ratio m of 0.2. 



Figure 2.- Vortex-blade intersection point path for various tip-speed ratios. Three-blade rotor; fiR = 


Tip- speed ratio 



1 I I 1 

10 280 320 360 




feet per second (198 meters per second). 












00 


O Flapwise bending mode 
□ Chordwise bending mode 



Figure 7.- Critical tip-speed ratios for various blade bending modes of cantilever blade. Three-blade rotor; QR = 650 feet per second U98 meters per second). 















M 

w 


Number of rotor rotations 

I- 1 -1- 2 -|- 3 -i— 4 -I 



Following-blade azimuth position, tp, deg 


Figure 12.- Blade-root (x = R) flapwise bending moment for single and repetitive Impulse. £ = 0; three-blade rotor; 
M ~ 0.3 ihilf = 1500); n = 3; Vj = 0.13; V 2 = 0.05; V 3 = 0.07. 
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